AD-A277  097 

liililHii 


•  1.  ag;>jc' ■  i  ■  s'.-i  1 .A' 

1  March  8,  1994  Reprint 

j  A.  sv^oTiTL-  TTnrTTTrri 

;  Current  Collection  and  Ambient  Plasma  Measurement  During 

Electron  Beam  Emission  at  Near  Geosynchronous  Altitudes  PE  62010F 

; _ _ _  PR  7601 

I  6,  AUTHOR'S)  '  TA  30 

!  ’  WU  06 

Shu  T.  Lai 


Phillips  Lab/WSSI 
29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010  .  ^  PL-TR-94-2046 


Reprinted  from  AIAA  94-0370  32nd  Aerospace  Sciences  Meeting  &  Exhibit 
January  10-13,  1994  Reno  NV 


Approved  for  public  release;  Distribution  unlimited 


This  paper  discusses  spacecraft  interactions  during  electron  beam  emissions  from 
a  spacecraft  rotating  in  sunlight.  We  present  an  improved  Langmuir  probe 
equation  for  the  current  collection  by  the  SCATHA  satellite.  The  equation  yields 
results  on  ambient  plasma  density  and  temperature.  The  results  compare  favorable 
with  the  very  few  measurements  obtained  by  other  methods  on  the  elusive  ambient 
plasma  at  near  geosynchronous  altitudes. 

i 

i 


14.  SUBJECT  TERMS 

Spacecraft  charging.  Spacecraft  interactions.  Electron  beam. 

Ion  beam.  Geosynchronous,  SCATHA 

15.  NUMBER  OF  PAGES 

6 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION  1  19.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE  1  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-0 1-280-SS00 


Standard  Form  298  (Rev  2-89) 

P'^cnoed  by  "ita  239* '8 
298-102 


AIAA  94-0370 

Current  Collection  and  Ambient  Plasma  Measurement 

during  Electron  Beam  Emission  at  Near 

Geosynchronous  Altitudes 

Shu  T.  Lai 

Phillips  Laboratory 

Hanscom  AFB,  MA 


S  f.  5 


32nd  Aerospace  Sciences 
Meeting  &  Exhibit 

January  10-13, 1994/  Reno,  NV 


For  permission  to  copy  or  repubiish,  contact  the  American  Institute  of  Aeronautics  Astronautics 

370  L'Enfant  Promenade,  S.W.,  Washington,  D.C.  20024 


4 


Current  Collection  and  Ambient  Plasma  Measurement 
during  Electron  Beam  Emission  at  Near  Geos3mchronous  Altitudes 
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ABSTRACT 

This  paper  discusses  spacecraft  interaaions 
during  electron  beam  emissions  from  a  spacecraft 
rotating  in  sunlight  We  present  an  improved 
Langmuir  probe  equation  for  the  current  collection 
by  the  SCATHA  satellite.  The  equation  yields 
results  on  ambient  plasma  density  and  temperature. 
The  results  compare  favorably  with  the  very  few 
measurements  obtained  by  other  methods  on  the 
elusive  ambient  plasma  at  near  geosynchronous 
altitudes. 

INTRODUCnON 

Spacecraft  potentials  float  with  respea  to  the 
ambient  plasma.  During  electron  beam  operations, 
a  spaceaaft’s  potential  varies  in  response  to  the 
emitted  current  At  equilibrium,  the  spaceaaft 
potential  and  the  ambient  current  collected  are 
related  by  Langmuir’s  plasma  probe  equation.  In 
the  laboratory,  the  current  collected  by  a  Langmuir 
probe  varies  according  to  the  applied  potential.  In 
space,  however,  the  current  is  the  driving  force  and 
the  probe  potential  is  the  response. 

In  ideal  geometries,  such  as  a  sphere  or  an 
infinitely  long  cylinder,  the  equation  is  well  known. 
The  geometry  for  most  spacecraft,  however,  is 
neither  a  sphere  nor  an  infinite  cylinder.  For 
improved  spacecraft  charging  modeling,  it  is  better 
to  determine  and  use  the  correa  Langmuir  probe 
equation  rather  than  taking  an  ideal  geometry 
equation  as  commonly  practiced. 

The  determination  of  the  Langmuir  probe 
equation  for  a  given  spacecraft  requires  data  of 
ambient  current  and  spacecraft  potential.  This  is 
often  difficult  because  the  current  data  often 
represent  the  sum  of  cunents  of  different  processes 
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and  the  potential  may  be  due  to  some  complex 
causes.  For  example,  when  a  spacecraft  is  rotating 
in  sunlight  [Lai,  et  al,  1986],  photoelectrons  leave 
the  satellite  surfaces.  When  elearon  beams  are 
emitted  from  a  spacecraft,  photoelearon  generated 
from  the  boom  surfaces  may  return  to  the  spacecraft 
body  [Lai,  et  al,  1987;  1989).  Magnetic  storms  can 
affect  spacecraft  potential  but  storms  fluctuate  and 
are  very  complex. 

In  this  paper,  we  determine  an  improved 
Langmuir  probe  equation  for  the  SCATHA  satellite 
by  uking  the  current  voltage  dau  taken  on  Day  70, 
1981  on  the  satellite.  SCATHA  was  launched  in 
1979  for  the  purpose  of  spacecraft  charging  studies. 
Its  shape  was  a  short  cylinder.  Its  length  was  about 
0.8  m  and  so  was  its  diameter.  Its  altitude  was  near 
geosynchronous.  It  was  equipped  with  two  50  m 
booms  (SCIO)  deployed  oppositely  in  the  equatorial 
plane.  The  "tip*,  farther  out  20  m,  of  each  boom 
was  coated  with  copper  beryllium  (CuBe)  [Aggson, 
St  aL,  1983].  The  booms  are  electrically  isolated 
from  the  satellite  body.  SCATHA  rotated 
perpendicular  to  sunlight  at  about  once  per  min. 


Figure  1.  The  SCATHA  satellite  and  its  SCIO  booms. 
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We  choose  Day  70  for  the  following  reasons. 
The  environment  was  quiet  During  the  period 
chosen,  SCATHA  emitted  an  electron  beam  at  a 
constant  energy,  300  eV,  with  the  beam  current, 
increasing  continuously.  As  beam  electrons  left,  the 
spacecraft  potential,  0^,  increased  and  therefore 
ambient  elearons  came  in  to  balance  the  cunent 
The  potential  was  governed  by  the  balance  of 
currents.  In  this  situation,  the  varying  beam  current 
was  the  known  driving  force  while  the  satellite 
potential  was  a  response. 

We  first  identify  the  various  physical  processes 
occurring  during  the  electron  emissions.  We  then 
use  the  data  of  the  Langmuir  probe  regime  to 
dete.^nine  an  improved  formula. 

CURRENT- VOLTAGE  DATA 

The  SCIO  potential  ^scio  ^  P''*®  ^ 

^scio  “  ^ 

which  is  the  difference  between  <(>^^  the  potential  of 
the  tip  of  a  boom,  and  0,  the  spacecraft  potential. 
The  boom  material  has  high  secondary  emission. 
Therefore,  when  the  environment  is  quiet  and  there 
is  no  severe  charging,  the  is  near  the  ambient 
plasma  potential  (i.e.  0^  »  o,  and  ^scio  ^  *  *^)' 

Figure  1  shows  the  satellite  and  its  booms  in 
sunlight.  Figure  2  shows  the  oscillating  0  data 
taken  during  the  electron  beam  emission  period. 
The  electron  beam  energy  was  300  V  constantly. 
The  satellite  was  rotating  perpendicular  to  the 
direction  of  sunlight  with  the  booms  in  the 
equatorial  plane.  The  booms  subtended  an  angle  $ 
with  sunlight  9  »  0°  or  180°  when  the  booms  were 
parallel  to  sunlight  The  dam  dropouts  were  due  to 
calibrations. 

The  maxima  of  the  amplitude  |  |  occurred 

when  the  booms  were  parallel  to  sunlight  (9  s  0° 
or  180°)  and  the  minima  at  =  90°  or  270°.  Both 
ihe  average  potential  (ft  and  the  amplitude  of 
oscillation  |  Max  (ft  •  Min  (ft  |  increased  as  the  beam 
current  increased. 


Figure  2.  The  SCIO  potential  and  beam  burrent  data. 

THEORETICAL  INTERPRETATIONS 

The  spacecraft  potential  (ft  is  determined  by  the 
balance  of  currents  according  to  KirchhoB*s  law.  In 
the  Langmuir  probe  model,  the  currents  are 
balanced  as  follows: 

for  e<^  >  kT.  In  eq(l)  le(<^=0)  and  lj((^=0)  are  the 
ambient  electron  and  ion  currents  collected  at  (ft=Q, 
T  the  ambient  plasma  temperature,  I,,  the  beam 
current,  I^  is  the  returning  beam  current,  I,  the 
secondary  electron  current,  Ip,,^  the  photoelectron 
current  from  the  spacecraft  body  and  Ip|,  the 
photoelecuon  current  from  the  booms  to  the 
spacecraft  body.  ^  =  1  for  a  spherical  spacecraft 
and  2/ir^^  for  an  infinite  cylindrical 

spacecraft  The  ambient  ion  current  Ij  is  usually 
two  orders  of  magnitude  smaller  than  1^  [Reagan,  et 
al.,  1980].  Therefore  Ii(<^)  is  negligible  unless  (ft 
reaches  thousands  of  Volts  negative.  In  eq(2),  (ft  is 
positive  and  increases  as  the  electron  beam  current 
I(,  increases. 
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We  contend  that  the  oscillatoiy  behavior  of  the 
potential  (F1g.2)  is  due  to  the  modulation  by  the 
photoelectron  current  Ip,,  from  the  booms  to  the 
satellite  body.  The  booms  were  near  ambient 
potential  while  the  satellite  body  was  at  a  positive 
potential  during  beam  emissions.  As  the  satellite 
potential  0  increased,  the  potential  sheath  extended 
to  engulf  part  of  the  booms.  Let  1.^  be  the  total 
photoelectron  current  generated  from  the  booms. 
The  fraction  of  Ip  going  away  to  the  ambient  plasma 
depends  on  the  extent  of  boom  engulfrnent  (Figure 
3).  When  «  =  0°  or  180°,  Ip^  =  0;  when  d  =  90° 
or  270°,  Ip,,  was  maximum. 

MATHEMATICAL  MODELING 

We  now  derive  Ip,,  and  /.  Figure  3  shows  the 
geometry  of  the  boom  angles.  The  angle  9  between 
the  sunlight  unit  veaor  =  (cosO,  sind,  0)  and  a 
boom  radius  unit  veaor  s  (0,  cos^,  sin^)  is 
given  by 

Sj.pj  =  cos  9  =  Sind  cos^  (3) 

For  a  surface  area  dA,  the  effective  sunlit  area 
is  dA  cos9,  where  dA  »  p  dr  d^  and  r  is  the 
distance  along  a  boom. 

Ipk  "  coa<p;^(0,(p)/[<|)(r)] 

(4) 


Figure  3.  Geometry  of  the  boom  angles. 


where  jph(d,v)  is  the  photoemissivity  of  CuBe,  the 
boom  surface  material,  and  L  the  length  of  a  boom. 
In  terms  of  surface  reflectivity  R  [Lai,  et  aL,  1986), 
written  as 

y^(0.»)  -4*(90*.0')  1 1  -  «(0.»)]  (5) 

Since  jp,,(d,9)  of  the  booms  has  not  been  measured, 
we  assume  a  constant  jp,,,  which  is  to  be  determined. 
If  both  jp,,  and /are  constants,  eq(4)  becomes 

l^~2frLJ^m9  (6) 


Figure  4.  Sheath  engulhnent  and  the  photoelectron  energy 
partition  function. 

We  now  calculate  /{<^(r)].  By  means  of  the 
partition  of  the  photoelectron  energy  distribution 
funaion  (Figure  4)  at  distance  r  along  a  boom,  the 
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(10) 


fraction  /  is  given  by 

rdEEexp(-ElkT^) 

-  (7) 

f’dEexpC-E/kT^) 

where  we  adopt  the  Debye  form  [WhippU,  et  al., 
1974]  for  the  satellite  sheath  potential  <^(r): 

^r)  -  (8) 

r*R, 

where  R,  is  the  radius  of  the  satellite  body  and 
0(rsO)  is  the  satellite  potential  0  (eq.2).  We 
integrate  eq(7)  to  obtain  a  handy  formula  of  the 
fraction  /  of  the  boom  photoelectron  current  going 
away  into  the  ambient  plasma: 

which  appears  in  the  integral  (eq.4).  In  eq(9),  when 
0  -•  0,  /  -*  fr,  when  0  -•  00,  /  -*  1  properly. 
Using  the  above  model,  Lai  [1994]  has  Otted  the 
data  of  Fimre  2  and  obtained  ^  3.S 
nanoamp/cm  ,  which  agrees  with  the  estimate  of 
Kellogg  [1980]. 

LANGMUIR  PROBE  FORMULA 

We  now  choose  the  potential  <f>  data  when  0  - 
0®  or  180®.  Since  Ip^  =  0  at  these  moments,  there 
is  no  photoelectron  current  involved. 

Three  regimes  can  be  identified  in  Figure  4.  In 
Regime  I,  the  spacecraft  potential  is  low. 
Therefore,  photoelecttons  and  secondary  electrons 
can  leave  the  spacecraft  body.  In  Regime  III,  the 
spacecraft  potential  is  near  the  beam  energy. 
Therefore,  beam  saturation  and  partial  beam  return 
can  occur. 

In  Regime  II,  there  is  no  photoelearon  or 
beam  return.  In  this  regime,  the  Langmuir  probe 
formula  (eq.2)  becomes 


where  a,  1^(0),  kT,  n  (between  1  and  1.1),  and  a 
(between  1  and  0.5)  are  to  be  determined.  In 
eq(lO),  when  ^*0,  we  have  1^(0)  »  I^.  Therefore, 
extrapolating  the  I-V  curve  to  0^0,  we  obtain  the 
ambient  plasma  current  1^(0)  »  10  /lA  during  the 
quiet  period.  This  1^(0)  result  agrees  with  the 
statistical  result  of  Purvis  et  al.  [1984], 
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Figures.  Solution  of  eq(ll). 

To  determine  the  exponent  a  (eqs  2,10)  and  the 
ambient  plasma  temperature  kT,  we  first  obtain  I^,’, 
the  differentiation  of  I,,  by  e<^  (Figure  S). 

Ii.l(,kT*e^)  (11) 

4  “ 

From  the  slope,  1/a,  and  the  intercept,  kVa.  we 
obtain  a  =  0.775  and  kT  *  23.2  eV. 

To  determine  the  ambient  plasma  density  n^, 
we  write  down  the  current  intercepted  by  the 
spacecraft  body  when  <f>  ^  0. 

^0)-ln^V,D‘l0]U  (12) 
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where  D  is  the  cross  sectional  area  of  the  satellite 
body,  the  electron  velocity  which  depends  on  IcT. 
Using  the  known  values  D  and  kT,  we  obuin  n^  s 
8.47  cm'^. 

If  value  of  lg(0)  (»  10  fiA)  is  accurate  and 
agrees  with  measurements  by  other  instruments,  we 
can  determine  n  by  uking  the  log  of  eq(lO). 

log/^-alog^l-»-^j>log^-»log//0)  (13) 

A  plot  of  logjo(Ij^  versus  logjQ(l+e<^/kT)  yields  an 
intercept  *  i  and  a  slope  *  0.775  in  Regime  II 
(Figure  6’).  Therefore  we  obtain  ^  I  and  a  » 
0.775. 


Figure  6.  Solution  of  eq(13). 


SUMMARY  AND  DISCUSSION 

By  modeling  the  various  processes  occurring 
during  a  period  of  elecuon  beam  emission  on  Day 
70,  1981,  we  have  identified  a  regime,  which  is 
properly  described  by  a  Langmuir  probe  model 
(eq.lO).  In  this  regime,  the  electron  beam  current 
controls  the  spacecraft  potential  without  beam 
return  or  photoelearons  involved.  We  have 
determined  the  parameters  in  the  Langmuir  probe 
equation  by  Htting  the  data  in  this  regime. 

The  results  obtained  are  (1)  the  boom  surface 
photoemissivity  jp„  »  3.5  nanoamp/cm^  (2)  the 
ambient  plasma  current  le(<^=0)  »  lO  fiA,  (3)  the 


Langmuir  probe  formula  exponent  a  »  0.775,  (4) 
the  Langmuir  probe  multiplicative  faaor  i,  (5) 
the  ambient  plasma  temperature  kT  «  23.2  eV,  and 
(6)  the  ambient  plasma  density  n^  *  8.47  cm’^. 

The  results  compare  favorably  with  the  very  few 
published  ones  obtained  by  other  methods.  The  j  ,, 
result  agrees  with  the  estimate  of  Kello^  [19%]. 
The  1^(0)  result  agrees  with  the  statistical  average  of 
Purvis  et  al.  [1984]. 

No  other  author  has  reported  a  method  of 
calculation  of  a  for  spacecraft  and  therefore  there  is 
no  comparison.  Our  a  result  (0.775)  is  between  1 
and  0.5,  corresponding  to  a  sphere  and  an  infinite 
cylinder  respectively.  Since  SCATHA  is  a  short 
cylinder  with  about  the  same  height  as  its  diameter, 
the  a  result  seems  reasonable. 

Again,  no  one  else  has  determined  fi  for  a 
spacecraft,  so  there  is  no  comparison.  Since  /x 
should  be  between  1  and  1.1,  corresponding  to  a 
sphere  and  an  infinite  cylinder  respectively,  our  /x 
value  (»1)  seems  reasonable. 

Whipple  [1981]  has  reponed  kT  «  64  eV  at  the 
geosynchronous  environment  during  a  quiet  period. 
Since  the  environment  is  unknown,  our  plasma 
temperature  kT  result  seems  reasonable.  Higel  and 
Lei  [1984]  used  ralaxation  sounding  to  obtain  an 
ambient  plasma  density  n^  »  3  to  7  cm'^  at  6.6  Rg, 
0500LT  and  2Kp  =  13.  SCATHA  was  at  7.3  Rg, 
05(X}LT,  and  =  14  -.  Our  ambient  plasma 
density  n^  (8.47  cm'^  is  comparable  on  the  high 
side. 
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